Particle-hole symmetric localization in optical lattices 
using time modulated random on-site potentials 
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The random hopping models exhibit many fascinating features, such as diverging localization 
length and density of states as energy approaches the bandcenter, due to its particle-hole symmetry. 
Nevertheless, such models are yet to be realized experimentally because the particle-hole symmetry 
is easily destroyed by diagonal disorder. Here we propose that a pure random hopping model can 
be effectively realized in ultracold atoms by modulating a disordered onsite potential in particular 
frequency ranges. This idea is motivated by the recent development of the phenomena called "dy- 
namical localization" or "coherent destruction of tunneling" . Investigating the application of this 
idea in one dimension, we find that if the oscillation frequency of the disorder potential is gradually 
increased from zero to infinity, one can tune a non-interacting system from an Anderson insulator 
to a random hopping model with diverging localization length at the band center, and eventually 
to a uniform-hopping tight-binding model. 



I. INTRODUCTION 

The interplay of quenched disorder and low dimension- 
ality in quantum systems has been a central theme in con- 
densed matter physics. A beautiful theory describes how 
the spectrum of disordered Hamiltonians follows from 
the symmetries of the problem 1 . In this work, we fo- 
cus on the one-dimensional non-interacting random hop- 
ping model, namely the tight-binding model with pure 
off-diagonal disorder. This model's properties have been 
extensively investigated theoretically using various tech- 
niques for many years. This model also exhibits physics 
manifest in several other models, such as quantum par- 
ticles connected by random strength springs, spin 1/2 
random XX chains, random quantum Ising chains in a 
transverse field, and random mass Dirac fermions. This 
model exhibits many surprising features. Early theo- 
retical work 2-5 focused on properties derivable from the 
mean local Green's function, notably the typical local- 
ization length which diverges as ~ \n\E\ and the mean 
density of states (DOS) also diverging as ~ l/\E(\nE 2 ) 3 \ 
as energy E approaches the band-center. Such behavior 
is very different from Anderson insulators where disorder 
appears in the diagonal terms of the Hamiltonian, and 
there are no singularities in the localization length or 
density of states spectrum. More recent work has stud- 
ied this model using real-space RG 6 and supersymme- 
try methods 7 and has uncovered more interesting results, 
most importanly, an additional length scale - mean local- 
ization length as a function of energy - which diverges as 
~ In |_E|. Most recently, the effect of random hopping 
amplitudes on interacting fcrmionic and bosonic systems 
has also been investigated. In the 2d fermionic case 8-10 , it 
has been shown that random hopping amplitude, on top 
of nontrivial spectral effects, could lead to a novel type of 
instability; in the bosonic case 11 ' 12 , a "Mott glass" phase 
has been predicted in addition to usual Mott insulating 
and supcrfluid phases. 



On the other hand, the experimental realization of 
a pure random hopping model has proven to be ex- 
tremely difficult. This is mainly because diagonal dis- 
order inevitably arises when we try to disorder bond val- 
ues, and any amount of diagonal disorder would break 
the particle- hole symmetry of the random- hopping model 
and thereby destroy the universal behavior of this class 
of disordered Hamiltonians. 

In this work, we propose that a pure random hopping 
model can be realized in optical lattices by first creat- 
ing an Anderson insulator and then modulating the dis- 
ordered on-site potential energy periodically. Our idea 
is closely connected to recent work on the phenomena 
dubbed "dynmaical localization" or "coherent destruc- 
tion of tunneling" in double wells, semiconductor super- 
lattices, and recently optical lattices 13-33 . While one can 
use this path to the random hopping model in any di- 
mensionality [d < 3), we will concentrate on its Id ap- 
plication below. 

The basic idea of dynamical localization is the follow- 
ing. Consider a particle in a double-well potential with 
a tunneling amplitude J and a time-modulated potential 
energy offset: 

H = -J(a t 6 + 6 t a) + Fcos(o;t)(a t a-6 t 6). (1) 

By performing a unitary transformation with 



U 



-%—■ sin(u;£)(a'a — b* b) 



(2) 



one can readily obtain that in the large-w limit, the orig- 
inal Schrodingcr equation is transformed to the effective 
Hamiltonian 



H, 



eff 



J Jo 



V 



(a^b + tfa), 



(3) 



where j7o is the zeroth Bessel function. Thus one can 
see that the effect of time-modulating the potential en- 
ergy is to simply renormalize the tunneling amplitude 
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FIG. 1: Phase diagram of the model (4) studied in this work. 
At zero frequency, the system is an Anderson insulator; when 
the frequency lo is comparable to the disorder width 7, the 
system behaves as a random hopping model; when w is much 
larger than 7, the system crosses over to the uniform-hopping 
tight-binding regime, which is fully achieved when ui = 00. 
We also present some interesting and puzzling results for the 
regime < uj <C 7. 



J in the large-w limit. This phenomena has been ob- 
served in experiments 27,28 ' 33 . It has also been proposed 
as a method to tune interacting bosons through the 
supcrfhiid-insulator transition 21,23 , to observe the ana- 
log of photon assisted tunneling and Shapiro steps 16,22 , 
and to manipulate the localization properties of Ander- 
son insulators 17,24 . Recent experimental work has con- 
firmed some of these proposals 34,35 . 



calization properties of the following model: 
H = H + 2Vcos(ut), 



N-l 



Ho = -J ^2 ( C 'l C ™+ 1 + C 1+1 C «): 



(4) 



iV 



v = 2> 



Cm C r . 



where N is the system size, and we assume the on- 
site potential v n obeys a uniform distribution between 
[—7/2,7/2]. In Sec. II, we will introduce the Floquct 
formalism with which we obtain the localization length 
and the density of states of the time-dependent Hamil- 
tonian (4). In Sec. Ill, an effective time-independent 
Hamiltonian is shown to emerge from the high-frequency 
limit of the original time-dependent Hamiltonian (4) in 
analogy to the dynamical localization phenomena. Nu- 
merical results from both the Floquet calculation and the 
effective Hamiltonian calculation are presented and dis- 
cussed in detail in Sec. IV. Next, we discuss possible 
experimental methods to modulate on-site disorder and 
to detect signatures of random hopping models in optical 
lattices. Finally, we summarize our results in Sec. VI. 



For our purpose, it suffices to notice that the orig- 
inal potential energy V resides in the renormalization 
factor of the hopping amplitude. Thus, if one modu- 
lates an Anderson insulator instead, one expects that 
the disorder in the onsite energy should be transformed 
into the disorder of hopping amplitude in the same way. 
In other words, one obtains the random hopping model 
by rapidly modulating the disordered potential energies 
of an Anderson insulator. As we will see in later sec- 
tions, however, if the modulating frequency uj is much 
larger than typical potential energy, this randomness in 
effective hopping amplitude is suppressed, and we obtain 
a uniform-hopping tight-binding model. Therefore, the 
random hopping model behavior survives when the fre- 
quency uj is comparable to the typical potential energy. 
In summary, as the oscillation frequency of the poten- 
tial energy is gradually increased from zero to infinity, 
one can tune a non-interacting system from an Anderson 
insulator to a random hopping model with diverging lo- 
calization length at the band center, and eventually to a 
uniform- hopping tight-binding model (see Fig. 1). 

Note, also, that a different but related model was stud- 
ied in Refs. (17,24,25), where an Anderson insulator with 
stationary disordered potentials but oscillating linear po- 
tential (e.g., uniform AC electric field) is considered. In 
that model, the localization properties of the Anderson 
insulator can be manipulated by the oscillating linear po- 
tential, but the random hopping model behavior is not 
accessible. 

In the remainder of this work, we will analyze the lo- 



ll. COMPUTATION OF THE DENSITY OF 
STATES AND THE LOCALIZATION LENGTH 

For a time-periodic system with Hamiltonian H(t) and 
period T = 2tt/uj, by the Floquet theorem, its wavefunc- 
tions can be written in the form 



m 



iEI 



<Kt), 



(5) 



where E is the quasienergy defined modulo w, and 
<j>(t + T) — </>(£). Here and in what follows we have 
set h = 1. This well-known result is the analog of the 
Bloch's theorem for particles in a periodic spatial poten- 
tial. To solve for E and <f)(t), one approach is to rewrite 
the Schrodingcr equation 



as 



id t1 p(t) = H^(t) 



H F (j) = E<j>, 



where Hp is the so-called Floquct Hamiltonian 
H F = H - idt 



(6) 



(7) 



(8) 



which is a matrix in the augmented Hilbcrt space T-L x 
T, where H is the original Hilbert space, and T is the 
frequency space 36,37 , and then to find the eigenvalues and 
cigenstates of Hp. Alternatively, it is also well-known 
that e~ tET and 4>(T) are the eigenvalue and eigenstate 
of the Floquet operator 



J- = T cxp 



dtH{t) 



(9) 



where 7" is the time-ordering operator. 

To obtain the density of states, we work with the lat- 
ter approach. First, we calculate the Floquet operator 
T by the numerical Trotterization procedure. Then, we 
diagonalize the Floquet operator T to find the quasiener- 
gies which we define to be in the "first Brillouin zone" 
—ui/2 < E < to/2. Then, we obtain the cumulative dis- 
tribution function of the quasienergies, average it over 
many realization of disorder, numerically differentiate it 
with respect to quasienergy, and finally obtain the den- 
sity of states. 

We would also like to obtain the localization length 
of this model for arbitrary frequency uj. For one- 
dimensional non-interacting time-independent systems 
with iV-sites, we recall that the localization length of a 
state with energy E is given by 3 



m = -^ M<h * m 



where the Green's function 



G(E) = (EI-Hy 



(10) 



(11) 



I is the identity matrix. 

Following Ref. 25,38, we generalize the concept of lo- 
calization length of a time-periodic system by defininng 
it as the localization length of the time-averaged wave- 
function. In terms of the Green's function, it is 



t4t = - I™ ±-ln\G 1N (E)\ 
X(E) n^oo N 



(12) 



where 



G=(Q = 0\G F {E)\n = 0) (13) 

and Gf is the so-called Floquet Green's function: 

G F (E) = (EI-H F )- 1 . (14) 

Here, if we denote the frequency operator Cl = id t and 
its eigenstates 



Cl\n) = nu)\n), 



(15) 



|fi = 0) introduced above is simply the eigenstate with 

71 = 0. 

Next, we discuss how to compute G = (0 = 
0\Gf(E)\Q = 0). For a Hamiltonian of the form 



from 



H = H + 2Vcos(ut), 



(EI - H F )G F = I, 



(16) 



(17) 



we insert the resolution of identity in the frequency space 
and obtain 



J2(m\(EI-H F )\p)(p\G F \n) = I{ 



(18) 



Thus, 

[(E + niUj)I — Hq] G mn — VGm+Ln — VG m -l,n = ISmn, 

where 

G mn = (m\G\n). (19) 

To solve Gqo(E) from this system of equations, we follow 
Ref. 25 to obtain 

G(E) = G Q0 (E) = (EI -H - V+ f - V- ff )-\ (20) 

where 



1 



eff E±lu-H - V i 



:V. 



TV V (21) 



E±2uj-H -V-V 

The number of iterations needed to ensure the conver- 
gence of V e f f is roughly proportional to 7/w. 

Finally, it would be interesting to mathematically 
check if the Thouless relation 3 between the density of 
states and the localization length holds for a Floquet sys- 
tem. 



III. EFFECTIVE HAMILTONIAN IN THE FAST 
OSCILLATION LIMIT 

In this section, we show that if the oscillation frequency 
lo is comparable or larger than the disorder width 7, the 
original Schrodingcr equation of a time-dependent Hamil- 
tonian can be transformed to that of a time-independent 
effective Hamiltonian. For the original Schrodinger equa- 
tion 

\dtip = Hip, H = H Q + 2Vcos(cjt), (22) 

we define 

i/} = Ui>, U = e - 2isin M) v /". (23) 

Then the Schrodinger equation becomes 

id t $ = H e ffii>, H eff = U f HU - 2Vcos(uit). 

Using the operator identity 



e ij J C n C i C g- ir / C !i C n 






(24) 



we have 



H eff = - J Y1 [( c « c "+i + c Lk c «) 

n 

V (-i) m J ( 2 ( Vn ~ Vn +^ 



cos(mujt) 



+ i(4 C «+l - 4 + l C n) 

x jr (-iy n+i j, (2 



sin(mu)t) 



For uj larger than or comparable to 7, the argument of 
the Bessel functions is comparable to or smaller than 1. 
Hence Jq dominates over other Bessel functions, and we 
obtain an effective time-independent Hamitonian 



Inverse Localization Length Inverse Density Of States 



H 



eff ~ - J 5Z ( C « C «+1 + c «+l c n) Jo 



2(v n - V n +l) 



(25) 
which is a tight-binding model with a site-dependent ef- 
fective hopping amplitude 



Jeff,n — J Jo 



2(v„ -v n+ i) 



(26) 



When uj is not too large, this is exactly the random 
hopping model we arc looking for, and it should exhibit 
behaviors such as diverging localization length and den- 
sity of states at the band center. We can compute the 
localization length of this effective Hamiltonian by using 
(10) and compare with the exact calculation using (10): 



1 



\{E) 

'eff( 



lim l-\n\G 1N (E)\ 

Af->oo TV 



,u(E) = (EI-H^ } )-\ 



(27) 



It is also straightforward to compute the density of states 
of this model numerically. 
However, when uj 3> 7, 



Jo 



2(v r , 



(28) 



regardless of the value of v n . Therefore in this limit 
the system behaves like a uniform-hopping tight-binding 
model (see FIG. 1). In this regime, we expect the local- 
ization length of every state to diverge, and the density 
of states diverges at the band edges instead. 



IV. NUMERICAL RESULTS 

We have computed the localization length and the 
density of states for various values of the frequency uj 
with fixed hopping amplitude J = 1 and disorder width 
7 = 10. The results from the Floquet technique for the 
original Hamiltonian (4) and those obtained from the ef- 
fective Hamiltonian (25) at uj = 3, 7, 35 are plotted in 
FIG. 2. One can see that when uj = 7 which is compa- 
rable to the disorder width 7 and when w = 35 which is 
much larger than 7, the results from the exact Floquet 
calculation and those from effective Hamiltonian calcula- 
tion agree quite well. At w = 35 3> 7 = 10, every state is 
completely delocalized, and the density of states diverges 
at the two band edges, as expected for a uniform- hopping 
tight-binding model. On the other hand, at uj = 7, both 
the localization length and the density of states diverge at 
the band center, which arc characteristic of the random 
hopping model, as expected. The case of uj = 3 is slightly 
more surprising: although the effective Hamiltonian docs 
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Inverse Localization Length Inverse Density Of States 

FIG. 2: Inverse of the localization length 1/A (a, c, e) and 
inverse of the density of states 1/p (b, d, f) vs. quasienergy 
E (in units of the hopping amplitude J). Disorder width 
7 = 10 J, hopping amlitude J — 1, oscillation frequency u = 
3 J, 7 J, 35 J. Solid lines are from the Floquet calculation of 
the original model (4); solid triangles are from the effective 
Hamiltonian H:J, (25). Averaged over 1000 realizations of 



disorder. 



l e// 



not work well, the system still exhibits diverging local- 
ization length and density of states at the band center. 
In FIG. 3, we plot the localization length for more values 
of uj from 3 to 35, and the trend from random hopping 
model behavior to uniform-hopping tight-binding model 
as uj increases is clearly seen. 



Near the band center, we fit the results of localization 
and the density of states to well-known analytical results 5 
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FIG. 3: Inverse of the localization length 1/A vs. quasienergy 
E (in units of the hopping J). Averaged over 1000 realizations 
of disorder. Hopping amplitude J — 1, disorde width 7 = 
10J, oscillation frequency uj — 3J, 5 J, 7 J, 9J, 15J, 35J. 



(see FIG. 4) 



p(E) = N 
\(E) = 



2a 2 



\E(ln(E/E f)zy 
2\ln(E/E' ) 2 \ 



(29) 



where N is the system size, E and E' are energy scales, 
and a is the standard deviation of the logarithm of the 
effective hopping amplitude square In J 2 et, with [see Eqn. 

(26)] 



J, 



eff,n 



JJo 



2(Vn 



(30) 



We can easily evaluate a numerically to be 1.535 given 
lu = 7, 7 = 10. Fitting numerical results of localization 
length and density of states, wc obtain 



<?fit,\ = 1-496, Ofit^p = 1.677, 



(31) 



which arc quite close to the theoretical value 1.535 ob- 
tained above, further confirming our expectation that 
random hopping model behavior can be achieved by fast- 
modulating the onsite energy of Anderson insulators. 

At frequencies much smaller than 7 and the original 
hopping strength J, interestingly, the system has quite 
large localization length in this regime. In FIG. 5, we 
plot the inverse of the localization length vs. the label 
(e.g., 1st, 2nd, 3rd, ..., 100th) of every Floquet eigen- 
state for uj = 0.01,0.05,0.1,0.5,1,3,7 with system size 
N = 100 (the total number of Floquet states equals the 
system size N) . One can see that from u> — 7, when uj is 
lowered, first the localization length decreases (inverse of 
the localization increases), but around uj = 3 this trend 
is reversed, and all the states become more and more de- 
localized at smaller frequencies. At w = 0.01,0.05 and 
0.1, all the states have almost equally large localization 
length. 
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FIG. 4: Fitting the inverse of the localization length A 1 and 
the inverse of the density of states p~ near the band cen- 
ter to their analytical form p(E) = N ■ 2a 2 /\E(ln(E/E ) 2 ) 3 \ 
and \(E) = 2\ln(E/E' ) 2 \/a 2 . Here, E a and E' are energy 
scales 5 , and the theoretical value of a is the standard devi- 
ation of \nj 2 ff, which is 1.535 for our choice of parameters 
here. The fitted value of a is 1.496 for 1/A and 1.677 for 1/p. 
Hopping amplitude J — 1, oscillation frequency ui — 7, and 
disorder width 7 = 10. 



This trend of derealization at small frequencies is 
quite puzzling, but it is interesting to notice that in a 
similar model where an Anderson insulator is manipu- 
lated with an AC electric field 17,24 ' 25 , an analogous de- 
localization trend was found. An intuitive argument was 
given in Ref. 25, where it is argued that the modulation 
with frequency uj allows electrons to absort or emit in- 
teger numbers of "photons" with energy uj. When uj is 
smaller, more states with quasienergy E ± nuj arc in the 
original energy band. Since the scale of the localization 
length of the new state should be set by the state with 
the largest localization length among all the states with 
energy E ± nuj, a smaller uj imples that it is more likely 
for the new state to have larger localization length. This 
intuitive picture could be of some relevance to our case 
here as well. 
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FIG. 5: (Color online.) Inverse of the localization length 1/A 
vs. the label (e.g., 1st, 2nd, 3rd, ..., 100th) of every Floquet 
eigenstate for u — 0.05, 0.1, 0.5, 1, 3, 7. System size N = 100, 
hopping amplitude J = 1, oscillation frequency uj = 7, and 
disorder width 7 = 10. 



EXPERIMENTAL FEASIBILITY 



two laser beams (one red and the other blue detuned) 
through the same speckle potential. 

The main experimental probe to detect Anderson lo- 
calization in cold atom system has been time-of-flight 
spectroscopy 40 ' 41 . In particular, for weak disorder when 
the condensate occupies a delocalized state, the con- 
densate ballistically expands when the confining poten- 
tial is removed. On the other hand for strong disorder 
potentials, the condensate occupies localized states and 
ceases to expand at a characteristic time after released 
from the trap. As we have seen, in the random hop- 
ping model some of the states are localized while oth- 
ers are delocalized (near the band center). Thus disen- 
tangling such behavior using time-of-flight spectroscopy 
alone might prove to be an experimental challenge. On 
the other hand, the well-developed technique of Bragg 
spectroscopy 42,43 allows the access to the spectral func- 
tion and therefore the density of states of quantum gas. 
Thus, perhaps the most promising way of detecting the 
Dyson delocalized state is through its distinct single 
single-particle density of states near the band center 
given by Eq. (29) measured through Bragg spectroscopy. 



The main motivation for our work is to experimen- 
tally realize the particle-hole symmetric disorder classes. 
Following the route outlined above, requires creating an 
Anderson insulator in optical lattices, and then modu- 
late the onsite disorder potential periodically in time. 
The behavior of ultracold atomic systems subject to dis- 
order, and the resulting localization phenomena is a field 
of intense current research 39 . Experiments have relied 
on two methods to introduce disorder into such sys- 
tems. The first involves using two incommensurate opti- 
cal lattice potentials, providing an effective realization of 
the Aubrey Andre model which has been shown to give 
localization 40 . The second method method uses a speckle 
potential produced by passing a laser through a diffusing 
plate which directly imprints a disorder potential 41 . 

The main challenge in realizing the phenomena intro- 
duced in this work is producing time-dependent disor- 
der potentials which periodically attract and repell the 
atoms in the optical lattice system. The most direct way 
to achieve this is to periodically change the detuning of 
the disorder potential potential from red to blue. That 
is, the disorder potential is given by 



V(r) 



3ttc 2 



t) 



(32) 



where c is the speed of light, ujo is the atomic resonance 
frequency, A = uj — ujq is the detuning frequency. Thus, 
the sign of the disorder potential can be periodically 
changed by periodically changing the detuning. This can 
be achieved by using an acousto-optic modulator to con- 
tinuously vary the laser frequency. However, sweeping 
through the rcsponance can produce undesirable atomic 
losses. Thus it might be best to periodically alternate 



VI. SUMMARY AND DISCUSSION 

The simple partical-hole symmetric class of quantum 
disorder problems exhibits many fascinating properties 
such as diverging localization length and diverging den- 
sity of states at the band center. It is difficult, however, 
to realize it experimentally because its crucial particle- 
hole symmetry is destroyed by any amount of diagonal 
disorder. Our work suggests a realistic path to effec- 
tively realizing random hopping models in optical lattices 
by fast-modulating the disordered potential energies of 
an Anderson insulator, but without need for fine tuning 
the local potential. Our idea follows naturally from the 
recently studied, so-called, "dynamical localization" or 
"coherent destruction of tunneling" phenomena as well 
as from recent efforts to observe Anderson insulators in 
optical lattices by using various ways to artificially pro- 
duce disordered potentials. 

The setup we suggest can explore more than just the 
static properties of a special disorder model. After all, 
we are describing the response of a quantum system to 
a strong time dependent potential. By exact diaonaliza- 
tion of the Floquet operators arising from our model, we 
explored the spectral properties of a Id system subject 
to strong time-periodic disorder. As shown in FIG. 2 and 
FIG. 3, the special features anticipated from the Id ran- 
dom hopping problem arise in a wide range of frequency 
modulations. Even with moderately small frequency uj 
where the effective Hamiltonian (averaged in the vector- 
potential-only gauge) does not provide a satisfactory de- 
scription, the localization length and the density of states 
still show the random-hopping-model behavior (see the 
case of uj = 3 in FIG. 2). Furthermore, a fit of our results 
near the band-center to the well-known theoretical form 5 



gives good agreement (see FIG. 4). 

Our model, however, gives results which we do not yet 
intuitively understand in the low frequency limit. As one 
can see from FIG. 3, the localization length gradually de- 
creases if the frequency uj is reduced from 35 to 3 given 
the disorder width 7 = 10. When the frequency u> is 
further lowered, our numerical results reveal that the lo- 
calization length starts to increase again, and it becomes 
quite large for w = 0.05 and u = 0.1 (see FIG. 5). We 
should recall that the terminus of the limit u — > is the 
well-known Anderson insulator, in which all states are 
localized. We intend to study the localization properties 
of such Floquet operators with strong disorder in future 
work. 

While here we concentrated on the realization of par- 
ticular classes of quantum disorder problems, the use of 
time modulated Hamiltonians could be a general path 
to the experimental realization of desirable systems. In- 
deed, this philosophy is already apparent in the quest 



to realize low-filling quantum Hall states in cold atom 
systems 44 , and we expect that in both cold-atom and 
solid-state based devices, time dependent Hamiltonians 
will become a standard tool for the stabilization of unique 
many body wave functions. 
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